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Abstract: We derive a generalised form of flow equations for extremal static and rotating 
non-BPS black holes in four-dimensional ungauged N = 2 supergravity coupled to vector 
multiplets. For particular charge vectors, we give stabilisation equations for the scalars, 
analogous to the BPS case, describing full known solutions. Based on this, we propose a 
generic ansatz for the stabilisation equations, which surprisingly includes ratios of harmonic 
functions. 
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1. Introduction 

The study of black holes in theories with eight or more supercharges resulting from string 
theory compactifications has proved to be a very useful tool in uncovering some of the 
structure of the underlying statistical systems. For supersymmetric black holes this task is 
facilitated by the fact that they exhibit the attractor mechanism and full supersymmetry 
enhancement near the event horizon |l|, |2|, ^. Using the constraints imposed by supersym- 
metry, general stationary asymptotically flat solutions have been found in ungauged N = 2 
Einstein-Maxwell supergravity, including higher-derivative corrections, both in four and 
five dimensions |Q, |5|, ^, ^, § • The spatial profile of scalars in these solutions follows a 
first-order gradient flow, which is integrable to (non-differential) stabilisation equations, ex- 
pressing the scalars in terms of harmonic functions. On the event horizon (the endpoint of 
the flow), the values of scalars are dictated by the charges through the attractor equations, 
independently of the asymptotic boundary conditions (the beginning of the flow).^ 

In contrast, when the requirement that the solutions must preserve some supersym- 
metry is abandoned, much less is known about the general structure of the supergravity 

^Some authors interchange the meaning of the terms "stabilisation equations" and "attractor equations". 
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solutions and the microscopic theory behind them. The simplest generalisation of BPS 
black holes to consider are the extremal black holes which do not preserve any supersym- 
metry (see fl^). These are known to share some desirable features with the BPS branch, 
most importantly the attractor phenomenon |11, 13 1. 

For theories with 8 supercharges coupled to vector multiplets in four and five dimen- 
sions,^ the general structure of these non-BPS extremal solutions is unclear, since only 
partial results are known. In the static case, a restricted set of examples can be found by 
simply changing the sign of a subset of the charges, which breaks supersymmetry |14, |l5|| . 
It was found that the non-BPS solutions exhibit flat directions in the scalar sector, in the 
sense that the scalars are not completely fixed at the horizon once the charges are chosen 
1 14]. However, these examples are not generic enough — they contain one less than the 
minimum number of parameters required for the most general solution to be derived from 
them by dualities. A solution that does contain enough parameters is called a seed solution. 

For cubic prepotentials, an appropriate seed was found in [16, 17] and the full duality 
orbit for the stu model was subsequently derived in ]18]. This full example clarifies how 
the non-supersymmetric solutions differ from their BPS counterparts in more than simply 
changing the signs of charges. In particular, they have flat directions that are subject to 
symmetries that act along the full flow, including the horizon ]19, 21]. 

If one allows for angular momentum, there are two types of single-centre extremal 
solutions which display attractor behaviour |2^. The over-rotating (or ergo) branch are 
very different from the BPS solutions, as they feature an ergoregion and are continuously 
connected to the Kerr solution ]p3| , 24, In contrast, the under-rotating (or ergo- free) 



black holes have a continuous limit to static charged black holes and seem to be tractable 
using BPS-inspired techniques. Recently, the single-centre under-rotating seed solution and 



various multi-centred generalisations were found in ]26, 27, 28, 30]. In these cases, the 



nontrivial parameter appearing in the static seed solutions can be viewed as the constant 
part of a harmonic function describing rotation. 

Despite the existence of these known solutions, finding an organising principle for their 
general structure has proven challenging. The best developed approaches are based on four- 
dimensional supergravity, where electric-magnetic duality limits the possible structures. 
One such framework is provided by the timelike dimensional reduction of Breitenlohner, 
Maison and Gibbons ]^l|, which relates black holes, regardless of supersymmetry (or even 
extremality) , to geodesies on the (augmented) scalar manifold. Given sufficient symmetry 
on the scalar manifold, solutions, including multi-centre black holes, may be generated 
with powerful group-theoretical methods, cf. ]32, 33, U, 35, |3^] and references therein. 
Unfortunately, this comes at the expense of the results being expressed less explicitly. 

A more direct perspective has been offered by the fake superpotential approach of 
Ceresole and Dall'Agata |3^]. They noticed that the rewriting of the effective black hole 
potential for the scalars |11] as a sum of squares is not unique, leading to more than one 
type of first-order flow for the scalar fields. The flow, which in the supersymmetric case is 
governed by the absolute value of the central charge, may be more generally controlled by 

^ Since the two are related by dimensional reduction, we do not make a distinction between them in this 
introduction. 
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a different function, called the fake superpotential. The derivation of first-order equations 
based on a superpotential has been subsequently extended to static non-extremal black 
holes and for a number of models superpotentials have been identified explicitly [p9| , 
41, 42, Q (see ||4^ for a synopsis of these developments and |45| for earlier related work). 



The superpotential method has been first applied to multi-centre black holes in ||46| , 
which directly generalised |^]. However, simplifying assumptions restricted the non-super- 
symmetric solutions, as in [33|, to threshold-bound configurations with mutually local 
charges and unconstrained relative positions of the centres. In view of the recent results 
on the integrability of the scalar equations of motion in black hole backgrounds [47, 48 1, 
one might expect also the more complicated multi-centre solutions mentioned earlier to be 
derivable from first-order fiows integrated to stabilisation equations. 

As a step in this direction, we study extremal under-rotating (ergo-free) black holes in 
compactifications of Type IIB string theory on Calabi-Yau manifolds, using the formalism 
of 1^] . In section ^ we relax the additional conditions of 1 46 ] to arrive at the general form 
of first-order flow equations for stationary extremal black holes. Unfortunately, we flnd 
that unlike their previously known special cases, to which they correctly reduce under the 
relevant assumptions, they generically do not lend themselves to explicit integration. 

In section therefore, we follow a bottom-up approach, trying to find stabilisation 
equations by rewriting known solutions (expressed in terms of physical scalars or affine 
coordinates) in a symplectically covariant way, using the projective (homogeneous) coordi- 
nates. We find that this is indeed possible for the under-rotating seed solution of ||2^, if one 
adds a ratio of harmonic functions to the standard vector of harmonic functions appearing 
in the stabilisation equations. Motivated by this, we introduce an ansatz for the general 
case that can incorporate all known extremal solutions. Our arguments are independent of 
the considerations in section |^, but the general form of the proposed ansatz is compatible 
with the generic first-order fiow equations. However, it is difficult to fully impose it in the 
rotating case. 

Finally, in section ^ we combine our general flow equations with the ansatz in the static 
case and connect to the fake superpotential formalism. Section ^ is devoted to concluding 
remarks, whereas in the Appendix we present a general heuristic argument justifying the 
presence of ratios of harmonic functions in the stabilisation equations. 



2. Bosonic action and special geometry 

In rewriting of the effective action as a sum of squares and deriving the fiow equations for 
stationary black holes in 4-dimensional N = 2 supergravity, we largely follow the method 
and the notational conventions of the two papers ||5|, ^ whose results we generalise (where 
we also refer the reader for more details and additional references). 

Omitting the hypermultiplets, which are immaterial for our discussion, the relevant 



bosonic action ^| 



54D = ^J^^{r*1- 2g^i dz" A *d# - A d) , (2.1) 



-3- 



contains neutral complex scalars (belonging to the vector multiplets) and abelian 
gauge fields (from both the gravity multiplet and the vector multiplets), all coupled to 
gravity. 

The scalars are affine coordinates on a special Kahler manifold, whose metric can 
be calculated from the Kahler potential K: g^i = dad^K, where da is shorthand for d/dz^. 

The field strengths are defined as = dA^ , where = (A°, A"-), a = \, . . . ,n^. The 
dual field strengths Gi are given in terms of the field strengths and the kinetic matrix Mij 

by 

Gi = luiNi.j i^F^ + ReNijF^ . (2.2) 

We will not need the explicit formulae in what follows, but both the Kahler potential 
and the kinetic matrix J\fij for the vector fields are derivable from a prepotential, F, which 
is a homogeneous function of degree 2. The prepotential itself is typically displayed in 
homogeneous projective coordinates (z" = jX^) and we will take it to be of the 
cubic type: 

F = —Dabc ^0 =■■ i^'ffi^) ' /(^) = -g^abc^" ■ (2.3) 

Surface integrals surrounding the sources of the field strengths and their duals define 
physical magnetic and electric charges, and qi, respectively: 

P' = ^ [ F', 11 = ^ I Gi- (2-4) 



47r J52 47r J52 

Prom a geometrical point of view, the above theory can be regarded as the bosonic 
massless sector of type IIB superstring theory in 10 dimensions compactified on a Calabi- 
Yau three- fold Mqy-^ The scalars of the vector multiplets parametrise the moduli space of 
complex structure deformations of Mqy- The complex dimension of this scalar manifold is 
given by one of the Hodge numbers of Mcy, = h^'^, with the Kahler potential K{z,z) 
being determined by the unique (up to rescaling), nowhere vanishing holomorphic (3,0)- 
form Ohoi, characterising Mcy- 



K = -\ii[i QhoiA^^hoi • (2.5) 

It will be more convenient later to work with the covariantly holomorphic version of the 
top form 

0(z,z) = e^(^'^")/2i^hoi, (2.6) 
whose Kahler covariant derivative reads 

Vn = {d + iQ)n, (2.7) 

where Q = lm{daKdz°') plays the role of the connection. In components: 

Van = dan + \daK n , Van = dan - \daK n = Q. (2.8) 



^Thanks to mirror symmetry, one can equivalently use the type IIA picture. 
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In the canonical syniplectic basis {a/,/?"^} for the third integral cohomology H^[Mcy,'^), 
we can expand as: 

n = X^ai - FjI3^ , (2.9) 

where the coefficients are the periods of the Calabi-Yau manifold with respect to the dual 
homology basis of three-cycles {A^ , Bj}: 

= [ n = [ n^p^, Fi= [ n = [ nAaj. (2.10) 

Jai JMcy JBi JMcy 

Fj are further identified with the derivatives of the prepotential F with respect to X^: 
Fj = dF/dX^ (we hope that no confusion with the gauge field strength two- form F^ 
arises) . 

Similarly, the five-form field strength T of the IIB theory, assumed to take values in 
ri^(M4) 0iJ^(McY,^), where 0?{M4) represents the space of two-forms on spacetime, can 
be written as 

T = F^ ®ai -Gi® 13^ . (2.11) 

By integrating the field strength over an appropriate two-sphere in space, we recover the 
charges as the coefficients of the thrcc-form V € H^{Mcy,'^)- 

V=^j^J = p^a:-q,p'. (2.12) 

The five-form F is self-dual in 10 dimensions, *io-^ = (*'X)o)J^ = F, where -k and o 
represent the Hodge operators in, respectively, spacetime and the internal CY manifold 
McY- A representation of the Hodge operator on the basis forms {a/, can be given in 
terms of a scalar-dependent matrix M.{J\f): 

{'f)^M-HU){^'], (2.13, 

\oaj) \ajj 

so that the selfduality constraint on F can be expressed in terms of components as 

Instead of the canonical basis, one may use the Dolbeault cohomology basis furnished 
by {Q,,Vafi,'DaO^,(},}, diagonalising the Hodge operator o on Mcy- 

on = -in, oVan = iVan, (2.15) 

and satisfying 

{n,n) = -i, {Van,Vj,n) = ig^i, (i?„q,q) = o (2.16) 

with respect to the antisymmetric intersection product 

{Ei,E2)= [ E1AE2. (2.17) 

JMnY 
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In this notation, the central charge Z can be written as 

z(r) = (r, 17) = qix' - p'Fi , (2.18) 

and, conversely, one can prove that 

(ri, Ta) = 2 Im[-Z(ri) Z{T2) + 5""^ VaZ{Ti) ViZ{T2)] , (2.19) 

where 

gab jg 

the inverse matrix of gj^. 
Another useful object is the symmetric Hodge product {Ei,oE2), which introduces a 
norm on H^{Mcy,^)- 

\\Ef = {E,oE). (2.20) 

An example of its utility in the context of the attractor mechanism is provided by the black 
hole (or effective) potential 

Kh = i||rf = -l(p q)M{M)(^^, (2.21) 

where we suppressed the indices on the charges. (In what follows we will often identify 
the elements of H^{Mcy) with the associated vectors built out of the components in a 
symplectic basis, such as T and {p^qj)"^ here.) The matrices M.(J\f) and M.{M) are 
related to each other by the symplectic metric X, 

M=1M, ^ = "o^) ' = "-^ ' (^-^^^ 

and are functions of the kinetic matrix Mij (see eg. [^]), although detailed expressions 
will not be needed in our considerations. 



3. Flow equations from the action 

In this section we derive generalised flow equations for non-BPS extremal black holes in 
Type IIB compactifications on Calabi-Yau manifolds, using the formalism of As these 
equations are not directly integrable, we describe an algorithm for solving them. 

3.1 The action as a sum of squares 

Since we are interested in asymptotically flat, stationary extremal black hole configurations, 
the ansatz that we use for the spacetime metric is: 

ds'^ = -e2^(dt + ojidx'f + e-2^5ijdxMx^' , (3.1) 

with the condition C/(x*), 0Ji{x^) — )■ as r = ^SijX'^x^ 00. The one-form lj = Wjdx* 
encodes the angular momentum of the system. 

The action is not invariant under electromagnetic duality rotations, but as re- 

marked in 1^, at the cost of discarding manifest Lorentz invariance, one can have a duality 
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invariant formalism [52, |5^. For this purpose it is convenient to introduce the following 
product of spatial 2-forms B,C G n^{R^) H^{Mcy): 



e^^ r 

(B, C) = ^ / B A [★o(oC) - AoC)w + *o{w A ★oC)] , (3.2) 

1 - JMcy 

where by we denote the Hodge dual with respect to the flat three-dimensional metric 
Sij and define w = e^^w. In general, boldface symbols will be reserved for quantities in the 
three spatial dimensions. 

With this notation, the effective action (p.lD can be written as 



54Dcfr = -T^ f dt f \2dUA-kodU-^e^^du}Airodu} + 2gTdz''Airodz^ + {J^,J^)]. (3.3) 
lovr J Jir3 L ^ J 

As shown in this action can be re-expressed as a sum of squares, giving first-order flow 
equations for stationary supersymmetric black holes, including multi-centre composites: 

^ - 2 Im ★oD(e-^e-'"17) + 2 Re D(e^e-'"17a;) = , (3.4) 
Q + da + ^e^^ *oduj = , (3.5) 



D = d + i(Q + da + ie^^^odw) , Q = Im (5„/^dz") . (3.6) 



where 



In the light of the considerations in |^8|, one might expect similar equations, involving 
a modified field strength jF in place of the actual field strength to exist for non- 
supersymmetric extremal black holes as well. In fS^ this was shown to be true for the 
special case when the fake field strength is related to the real field strength by a constant 
symplectic matrix. Building on this, we have found a more general way of writing the 
action as a sum of squares resulting in non-BPS first-order fiow equations, based on a 
non-closed and a new auxiliary one-form rj related to the non-closure of JF^. These two 
objects are constrained by two new equations which need to be satisfied in addition to the 
flow equations we obtained. 



While they will be derived and explained below, for ease of comparison with (3^, p.5| ), 
we state our non-BPS equations now. The first two equations are very similar to the BPS 
ones: 

:F - 2 Im ★oD(e-^e-'"0) + 2 Re D(e^e-'"J7a;) = , (3.7) 
h 



Q + da + r] + ie^^^odw = , (3.8) 



with ^ replacing ^ in the first equation and rj shifting the second.^ In addition, the two 
equations constraining our two auxiliary variables are: 

(J^, J^) = (:f, T)+2d'nAw, (3.9) 
77 A Im(^, e^e-^^r?) = (d^, Re(e^e-^°r?)) -\dr]Aw, (3.10) 



D and Q are defined as in (3.6) 
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where 

G = T-2lm ★oD(e"^e"'°0) + 2 Re D(e^e"'°Oa;) . (3.11) 

Now onto the derivation. In Q the crucial step in obtaining the first-order manifestly 
duality-invariant flow equations, solving second-order equations of motion, is to appropri- 
ately pair the derivatives of the scalars with the gauge fields and use the invariant product 
(|3.2| ) to re-express the Lagrangian. It was found in ^ that a good choice is 

g = T-2lm ★oD(e-^e-^"Jl) + 2 Re D(e^e~'"Jla;) . (3.12) 



Similarly to ji^], we generalise the above to eq. ( 3.11 ) by replacing the actual field 
strength with a two-form valued 'fake' field strength, ^ G r2^(]R^) (g) H^{Mqy), but here 
we define it by demanding only that it reproduces the original electromagnetic part of the 
action 

{T,T) = {P,T)-S, (3.13) 

up to a possible extra term described by the three- form S. The form of S will be de- 
termined at the end of this subsection by consistency arguments. Unlike the real field 
strength, we do not require P to be closed. 

One can then rewrite the Lagrangian in terms of Q as 



C = {G ,g) - A{Cl + da + -q + ie^^^odw) A lm{G, 
+ d [ 2'u; A (Q + da) + 4 Re(^, e^e"'"!!) ] , 



(3.14) 



provided that the new one-form 77, which needs to be introduced due to the possible non- 
closure of iP, satisfies 

r? A Im(a, e^e-^"0) = (dJF, Re(e^e-'"1])) + \S . (3.15) 

Because adding a total derivative to the Lagrangian does not change the equations of 
motion, one finds that equation ( 3.15| ) needs to hold only up to a total derivative. Observe 
that the phase a = a(x) is a priori an arbitrary function. 

A sufficient condition for a stationary point of the action (hence, for the equations of 
motion to be satisfied) is met by neglecting the boundary terms and requiring that the 
variations of the first two terms in ( 3.14| ) vanish separately, which leads to (J^^, |3.8| ). From 
(|3.8|) we obtain D = d — ir/, which substituted into (3.7) gives: 



:F-2Im[*o(d-ir7)(e-^e-^"0)] + 2Re[(d - ir7)(e^e-^"0w)] =0. (3.16) 
Differentiating one finds: 

d*od Im J7 - d (★o^? ReClj - d (r] A wlmClj = \dP , (3.17) 

where (l = e-^e-'^O.^ We see in particular that, as mentioned earlier, the fake field 
strength is not necessarily closed. 



^Note that in |5| the hat symboHsed what we call o here. 
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It is now possible to derive S in terms of the other quantities appearing in the rewriting. 
The fundamental observation is that in the dynamical system we are considering, the 
electromagnetic part of the Lagrangian acts as a potential for the remaining fields and, 
for solutions of the equations of motion with vanishing action, it is expected to equal the 
kinetic energy. So, by expressing ^ through equation ( |3.16| ) , one can compute (^, ^) and 
find: 



2dU A ★odf/ - ie^^dw A ★odcj + 2g^idz'' A ★odz" 



+ e^^du; A ★oda; + 2dQ A w 
= {J^, + e^^dw A ★odw + 2dQ A w , 

from which, using ( 3.13| ), (|3.8D and integration by parts, one obtains 



(3.18) 



-2d?7 A w . 



(3.19) 



Finally, substituting ( PI) in ( PI ) and ( p5|) leads to (U, |3lc| ). 

In summary, we have obtained a non-supersymmetric generalisation ( |3.7| , |3.8| ) of first- 
order equations . The generalised equations are expressed in terms of a fake field strength 
constrained by ( ^.13 ) to reproduce the original gauge part of the action. The non-closure 



of J- necessitates the introduction of a new, compensating object, rj, in eq. ( 3.15 ). The 
auxiliary three-form S appearing in (|3l^ ) and ([SlsD can be expressed in terms of other 
quantities through eq. ( p. 19]) . In comparison with the supersymmetric case we thus have 
two more unknowns, IF and t], constrained by ( p.9| , 3.101 ). Since they are mutually related, 
in any model for which F can be obtained by other means (as in section |5.2| ), r] can be 
eliminated as well. 

3.2 Solving the equations 

Whenever dF = and rj = 0, equation (|3.17|) reduces to the Laplace equation 

2d*odImA = 0, (3.20) 

which can be integrated to so-called stabilisation equations: 

2lm{e-^e-'"n) = n. (3.21) 

These express the period vector in terms of (possibly multi-centred) harmonic functions Ti 
throughout the flow. 

Using this result back in (3.16) one finds that the fake field strength is given by 



irodn-2d{e^^Renoj). 



(3.22) 



This is manifestly true for the supersymmetric case |Q, ^ and for the non-BPS setting 
considered in [46|, for electric {p^,0;0,qa) or magnetic {0,p'^; qo,0) charge configuration 
and axions Rez" set to zero.^ In the BPS case, one has J- = J- hy assumption, so that 



By charge redefinitions one can generate physically equivalent solutions also for other charge configu- 
rations |b4|. 
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the vector of electric and magnetic charges is determined through ( |2.12| ) to be 

r = -^ / irodn, (3.23) 

47r Js^ 

or in other words equal to the poles of the harmonic functions H. For non-supersymmetric 
solutions, J- is related to ^ by charge sign reversals, and the same holds for the poles of 
the harmonic functions T-L in ( 3.22| ) compared to the physical charges. 



In its general form, however, equation ( 3.17 ) cannot be solved directly, since the period 
vector (scalars), rj, uj and J- are all unknown and constrained by ( p. 10 ). A way out of 
this problem is to first make an ansatz for ImO, try to solve it for and U, find t] and uj 
from (|3.8| ) and then, by using ( 3.16| ) as a definition for IF, check if (3.9, 3.10) are satisfied. 

Let us see explicitly how to do that: we start by making an ansatz of the type 



2Im(e' 



(3.24) 



with J' a vector containing all the parameters in terms of which the solutions will be 
expressed. The solution for the components of (and hence for scalars) can then be 
obtained in the same way as the solutions to supersymmetric stabilisation equations ( |^ , 
see also section 2). We will indicate all quantities calculated with the aid of the ansatz 
by adding the subscript J'. 



We then proceed differentiating both sides of ( 3.24 ) and subsequently taking the inter- 
section product with the real and imaginary part of and with Va^^. With the definitions 
^ := -{dj,n) and := -{dJ,Van) one obtains: 



dU 
da + Q 
dz' 



Re(e-^°^') 
e- Imfe-'"^') = 



-e 

JJ - 



IJ2U 



(3.25) 
(3.26) 
(3.27) 



Note that ( |3.25| ) and ( |3.27 ) are the fiow equations for the warp factor and the scalars, 
while (|]2|) gives an explicit relation between a and the other quantities appearing in the 
rewriting. More specifically, ( 3.26| ) combined with ( |3.8| ) eliminates a, giving: 



{dJ,J) = 2fi-^^ T] + -koduj . 



(3.28) 



If we make an ansatz also for the angular momentum of the black hole cjbh (which must 
be expressed in terms of parameters appearing in J), we arrive at an expression for rj: 



2^ 



({dJ,J) - ★od^bh 



(3.29) 



Independently of the ansatz, we can use ( p. 16 ) as a definition of J- and then substitute 
(3.17) and (3.7) in (3.10). The left-hand side is clearly zero, whereas on the right-hand 
side we have an intersection product that we know how to compute. Neglecting the total 
derivative results in: 



(d^, Re(e2^17)) = (2 d*od Im VL, Ke{e^^ Vl)) - A ★or? - e^^r? ^duJ + \r] ^dw . (3.30) 
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The last term here cancels the last term of ( 3.1C1| ). This means that once all the variables 



have been expressed in terms of the parameters in the vector J ^ the consistency of the 
ansatz with the first-order equations ( |3.7| , and the constraint ( |3.1C1| ) can be verified by 
checking whether the following equation is satisfied: 

e2^-7 (d*od:r, Re ^j) = r,j^^ A ★or/j,^ + e^^^Vj^^ A dwbh • (3.31) 

As this is an equation for J, in principle it determines an ansatz that satisfies (^.7]), 



and ( 3.10 ), even though in practice one would not solve it for regarding it instead as 
a check for the specific form of an ansatz assumed beforehand. One would also still have 
to ensure that (3.9) is satisfied, which can be a rather non-trivial task. Finally, not all the 



parameters in the ansatz may be constrained by the equations of motion but should rather 
be fixed by appropriate boundary conditions. 

In the next section we will discuss known black hole solutions which satisfy ( 3.31| ) with 
a nontrivial and propose a generic ansatz. 

4. Stabilisation equations from an ansatz 

An important result for BPS black holes is the direct integrability of the first order fiow 
equations to stabilisation equations, even for multiple centres. As described in the be- 



ginning of section 3^, this result can be extended to some non-supersymmetric solutions. 
These examples, however, are not generic, in the sense that applying dualities on them 
does not lead to the most general non-BPS solution. 



Given the fiow equations in section 3.1, one expects to find a nontrivial rj and ^ in 
the general case. The non-closure of J-" implies that the corresponding expression for the 
period vector $7 in terms of charges and integration constants should be an anharmonic 
extension of ( |3.21| ), which must still be consistent with symplectic reparametrisations. To 
gain intuition about the possible terms, one can follow a bottom-up approach. Therefore, 
we consider known explicit solutions and rewrite the physical scalars and the metric in 
terms of Vt, aiming towards a generic ansatz that covers all single-centre solutions. 

4.1 Special solutions 

In order to be as general as possible and to minimise the ambiguity introduced in the pro- 
cess, we find it illuminating to start with a rotating black hole solution, so that the presence 
of an extra harmonic function describing angular momentum can provide guidance. Con- 
sider the rotating extremal black hole of [p^], which can be used as the seed solution for 



four-dimensional under-rotating black holes |23, 24, 25 1 in theories with cubic prepoten- 



tials. This is an almost BPS [^] solution in five-dimensional supergravity described by the 
harmonic functions: 

H^ = h^ + ^-, Ha = ha + ^, M = h+i^, (4.1) 

ly* ly* ly* ^ 

where , ha, b are constants that are related to asymptotic moduli, —p^ is the Kaluza- 
Klein magnetic charge, qa are the electric charges and J is the angular momentum of 



- 11 - 



the solution. Therefore, the associated four-dimensional charge vector is defined by the 
harmonic functions: 

Hc = [-H\^-Q,Ha), (4.2) 

whereas M controls the angular momentum and is invariant under symplectic transforma- 
tions. 



Using the 4D/5D dictionary of |57, 58, 56 1, one can rewrite the full solution given 



in five-dimensional notation |^7[] in terms of variables natural from the four-dimensional 
perspective. The metric is as in while the resulting expressions for the gauge fields 

and scalars in our notation are: 

J^ = *od?^c-2d(e2^ReOw), 2 Im 1^ = J = H + 7e , (4.3) 

where we use again the shorthand $7 = e~'^e~'"0 as in section |^. J is written in terms 
of a harmonic part, and a part containing ratios of harmonic functions 7^, which are 
respectively given by: 

= (i/O, 0; 0, if„), 7e= (O, 0; O) . (4.4) 

Finally, the metric functions are given by: 

■koduj = dM, e-2^ = i(f2, (l) = ^h{n) - M2 . (4.5) 

Here, /4 is the quartic invariant that appears in the entropy formula for cubic prepotentials 
(see 1 59 1 for explicit expressions) and the physical scalars are given by = X"" /X^, as 



usual. 

The expression for the gauge fields in ( |4.3| ) parallels the form of the BPS solutions 
( 3. 221 ) , differing in that the vector of harmonic functions associated to the physical charges 



|4.2| ), is related to the one appearing in the scalars by a single sign fiip, similar to [15|]. The 



period vector Vt is again determined through stabilisation equations similar to ( 3.21 ), so that 
the scalars are given in terms of harmonic functions describing the fiow from asymptotic 
infinity to the horizon. In particular, the asymptotic values of the scalars are controlled by 
the constant parts of the harmonic functions % and M, whereas the attractor equations, 
obtained in the limit r ^ 0, are controlled by the charges and the angular momentum p2| . 

The novel addition to is a ratio of harmonic functions that was not present in 
previous attempts to write non-BPS stabilisation equations and allows for nontrivial axions. 
Deferring the comparison to the rewriting of section |^ for the next section, we note that 
this solution leads through ( p. 29 ) to a nontrivial 77, given by 



.7 = e2^(d7^,^) = -e2^i^°d(^), (4.6) 

which demonstrates how ry is related to the anharmonic part of the solution. 

In the static limit, it is possible to show that all constraints and fiow equations of the 
previous section are indeed satisfied, if J in ( |3.24| ) is identified with the one in ( [4.31) . The 
presence of a nontrivial rj that follows from the ratio M/H^ through ( [4.6[ ) is crucial in 
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this respect. In the rotating case, we have verified the constraint ( 3.31| ), but it is more 



chahenging to verify the flow equations and especially the first constraint ( |3.9| ). 

As this is the seed solution for under-rotating extremal black holes, one can apply 
duality rotations on the full rotating solution using the stabilisation equations ( |4.3| ) to find 
the most general solution. Imposing that the angular momentum harmonic function M is 
invariant under duality transformations, the result is that a ratio of harmonic functions is 
generated in all other cases as well. 

For example, in the case of the stu model, one can explicitly dualise to the frame with 
only two charges present, corresponding to a D0-D6 brane system in Type IIA theory. For 
this model, the prepotential is as in ( |2.3[ ) with Dabc = \£abc\ and the scalar sector is then 
described by the choice (no sum on a = 1, 2, 3): 

^ = i.if'^ ; Ho, A'^if.) , 7^ = i^(l, ; -1, -A«) , (4.7) 



where 



Hi = h, + '^, H^ = \{ho + Y.ha) + '^ (4.8) 

a 

= -X^Hj , = -X^H+ , DabcX^X'^y = A^ , (4.9) 

= h{n) - m2 = {HoH^f-M^, (4.10) 

and A^ must be a constant. Note that the individual constants A'* appear only as multi- 
plicative factors in Ti and TZ, but not in e"'^, which depends only on the physical harmonic 
functions Hq and H^. It follows that the metric and gauge fields depend only on the com- 
bination A^, so that two of the A" correspond to fiat directions. The structure in (|4.7|) is 



consistent with the results on D0-D6 attractors in [S0| and seems to be generic for D0-D6 
solutions for all cubic prepotentials. 

It is interesting to note that unlike in ([4.3|), the harmonic part H is not related to the 
charges by sign flips, as one might expect. In fact the electric solution is special, in the 
sense that the flat directions can be described through (|4.3| ) by simply allowing for the 
missing harmonic functions to be constants, at the cost of making TZ more complicated, 
but still proportional to a single ratio as in the D0-D6 case. On the other hand, for both 
solutions the angular momentum harmonic function can be invariantly characterised by 
M = {T-LjTZ). The flat directions described by the A" are zero modes of this equation. 



4.2 The ansatz 

On account of the above observations, it is natural to propose an ansatz for the period 
vector that contains harmonic functions and ratios of harmonic functions, disregarding the 
precise relation to the physical charges, which is to be fixed later. In fact, it is simple to see 
that imposing consistency of any generic ansatz Im Q ^ H + TZ, leads to inverse harmonic 
functions. Since one can compute {7i,TZ) in two ways: 

{n,TZ) = lm{Cl,TZ) = -lm{Cl,n) , (4.11) 
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where H and TZ are a priori independent, it follows that {Ti, TZ) must be a scalar-independent 
quantity. The only other fields in the system are the scale factor and the rotation form cj in 
the metric, but since {H, TZ) does not carry a scale, ^ it cannot depend on e^, in accord with 
the explicit solution above, where {T-L,TZ) = M. In the static limit M reduces to a constant, 
in which case the constraint could be solved even if TZ were harmonic, but in the rotating 
case one has to reproduce the full function M, which depends on the angular coordinates. 
This implies a structure as in ( [4. 3D , with the anharmonic part, which then must be present 
even when the angular momentum is turned off. A more extensive argument about the 
kind of zero modes allowed for the scalars, which leads to the same conclusion, is given in 
the Appendix. 

Based on the linearity of symplectic reparametrisations and the fact that (^^) and ( [4.7| ) 
are seed solutions, we expect the structure seen in the previous section to be universal for 
all under-rotating extremal black holes. In other words, we take the point of view that 
there is no essential difference between static non-supersymmetric and under-rotating black 
holes, since they are continuously connected by setting to zero the nonconstant part of a 
single harmonic function, as in ( |4.1| ). Therefore, we propose the following form for the 
stabilisation equations for the scalars and the angular momentum: 

2Imf] = 2Im(e"^e-'"0) = + 7e, (4.12) 
i^odoj = {dTi,'H) + d{Ti,TZ) , (4.13) 

where Ti is & vector of harmonic functions and 7^ is a vector of ratios of harmonic functions. 
The integrability condition of the last equation implies that their symplectic inner product 
{T-L,TZ) is a harmonic function, while the scale function of the metric is given by: 



e-^'^ =i{n,n) = ^^h{7i + TZ) . (4.14) 

Note that when 7^ = and the charges carried by H are identified with the physical 
charges, one recovers the BPS stabilisation equations, as required. More generally, for a 



physically reasonable solution the harmonic and inverse harmonic functions in (4.12) are 
quite restricted due to various consistency constraints, both generic and based on known 
explicit solutions. The rest of this section is devoted to a discussion of these generic 
constraints and some of their implications. 

A first requirement is that in the near-horizon limit the scale factor e~^^ of an under- 
rotating black hole must reduce to [p^j: 



e 



oc |/4(r)| - J^cos^e, (4.15) 



where /4(r) is the quartic invariant of the model and J is the angular momentum. In 
the simple case of vanishing angular momentum, TZ is proportional to inverse harmonic 
functions and thus vanishes near the horizon. Therefore, a harmonic piece must always 



be present in the right hand side of (4.12), to make sense of the static solution in the 



^Here we refer to the symmetry of (3.14) under e*^ — ^ e^e^ , Vl e^Q,, gtj — > e^^ gij for constant D, 



inherited from the full conformal formulation of the theory l6l| 
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near-horizon region. Similar comments then apply for the full rotating case, hence it is 
impossible to have a physical solution for the scalars based purely on inverse harmonic 
functions. 

Going over to the constraints posed by the form of the full solution, observe that in 
the (necessarily static) BPS case the full scale function is simply e~^^ = I/^iT-L), where 
the charges are replaced by their corresponding harmonic functions. Similarly, for the stu 
model, where the most general non-BPS static black hole was explicitly constructed in |18| 



using the seed solution of |16, 17 1, it has been shown that the scale factor is shifted as 
e~^^ ~ Ii{'H) — where 6 is a constant that does not depend on the charges. 

Interestingly, for the known under-rotating seed solution the expression for in 
( |4.5| ) can again be found from ( 4.15| ) by replacing the charges and angular momentum by 



harmonic functions. Moreover, the additional constant b of 17, |18|] is identified with the 
constant piece in the harmonic function for the angular momentum in ( [4.5D , as in |2(:, 27 1. 
Therefore it is reasonable to expect that generically the scale factor is a function of the 
harmonic functions for the charges and angular momentum, thus allowing for the presence 
of a possible residual constant in the static solutions, when J is set to zero. 

Now, for an ansatz of the type ( [4.12 ) to describe the known solutions, the vector TZ 



must be such that ( 4.14 ) is consistent with the above comments, in particular with ([4.5|), 
so that 

= h{n + iz) = h{n) - {n, nf . (4.16) 

This equality poses very strong restrictions on TZ, as it does not appear in linear, cubic or 
quartic terms. In particular, the components of TZ must be such that Ia{TZ) and its first 
derivatives vanish, implying that it must have at most as many independent components 
as a two-charge small black hole. Then, given % and a model in which /4 is known, the 
linear term in IZ should vanish, further restricting its independent components. Indeed, IZ 
appears to have only one independent component in the explicit solutions (4^) and ([4.7|). 



For symmetric cubic models this can be made more precise, by Taylor expanding the 
left hand side of ( [4.16|) explicitly. In these models, the quartic invariant can be rewritten 
in terms of the central charge as 



where 



(ii -i2)^ + 4u-«5, (4.17) 



h = ZZ, 12 = g'^'VaZV-.Z , (4.18) 

i3 = lRe{ZN3iZ)), u = -llm{ZN3{Z)), (4.19) 

i5 = 9''^D,hcD-,i,g'''g'''g''''g'''D^Z VaZV^Z V,Z , (4.20) 



and 



iV3(Z) = D^^.g^'^g^^'g'^^VaZ V^Z %Z . (4.21) 
Although these five invariants all depend on the scalar fields and the charges, the combi- 



nation in ( 4.17 ) is scalar independent. In this form, it is easy to expand I4{T-L + TZ) and 



separately consider the different terms, since Z and its derivatives are linear in the charges. 
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Furthermore, as shown in ||62[, there are relations between the invariants above when the 
charge vector is that of a smaU black hole. The previous discussion suggests that TZ should 
have only one independent component, so we assume that it lies in a doubly critical orbit, 
in which case 

i2(7e) = 3^l(7^); i3(7^)=0; u{n) = 2il{n) ; W = 12i?(7^) . (4.22) 

A straightforward expansion of ( 4.17] ), using (4.22) leads to 

hin + n) = h{n) + {5h{u),n) - {n, nf , (4.23) 

where 6l4{7i) denotes the derivative of hiTi) and the identity ( 2.19| ) was used. Thus, the 
quadratic term reorganizes itself in the desired form without further assumptions.^ For a 
given model, TZ can then be determined by demanding that the linear term vanishes. 

This requirement is enough to ensure that the ansatz ( 4.12| ), together with the above 
assumptions, automatically satisfies the constraint ( 3.10 ), as we now show. First, note that 
for the ansatz in ( 4.121) , rj takes the form 



T] = e"^^ {dn,n) 



(4.24) 



as in (|4.6| ). In section 3^ it was shown that the constraint ( |3.1[l| ) is equivalent to (|3.3lJ ), 
which in view of the last result reads 



(d^odj^. Re n) = e^^ (d7^, H) A (★od?^, TZ) . 



Using ( 4.14 ) and ( 4.16| ), one can then show that 
(d*odJ,Rel^) = ^e^^{TZ,n) (d★od7^,?/) 
where we used the identity |]65|| 

Refj = 



{TZ,n) {irodTZ,dn) , 




-2U 



(4.25) 



(4.26) 



(4.27) 



and the last step follows from the fact that {TZ,7i) is a harmonic function. Finally, since 
TZ depends only on a single ratio of the form {T-L,TZ)/H (with H a harmonic function, 
cf. (4^)), it is possible to show that 

{n,TZ)dTZ = -{dTZ,n)TZ. (4.28) 

Combining the last relation with ( 4.26| ), the constraint (4.25) is identically satisfied, so 
that ( [4.12 , 4.13| ) is a solution of the constraint ( |3.1[1| ). The seed solution ( |4.3| ) also satisfies 
these relations by construction. 

This is a rather nontrivial result, as (3.31), due to ( p. 29 ), is a quartic equation for J'. 
Assuming this to be the general solution, the only constraint remaining at this stage is 



* Conversely, the decomposition in Im] can be used to show that the quartic invariant takes this form 
only if TZ lies in a doubly critical orbit [M. 
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( |3.9| ), which generahses the constraint on the fake superpotential for static black holes ||38|] 
to the case of under-rotating black holes. However, it is difficult to verify (|3.9[ ) and ( 3.16| ) 



explicitly for the seed solution (4.3), or find the general solution. In the next section, we 



give a more detailed comparison to the ansatz (4.12, 4.13D in the static limit. 



It follows that the only object missing for a complete characterisation of the ansatz for 
extremal solutions is an explicit form for % given a vector of physical charges. In view of 
the flow equations in the previous section, that would be equivalent to solving (|3.9[) which. 



upon using ([4.12| , |4.13D to determine the scalars and 77, becomes a quadratic equation for 
the physical charges in terms of Ti and TZ. 

Unfortunately, solving this constraint is not a straightforward task. The only a priori 
requirement on % is that it must be "BPS" in the sense that Ii{'H) >0 and that its quartic 
invariant should be related to the one of the physical charges by a sign flip. In fact the result 
should not be unique, as one can expect in view of the non-uniqueness in the rewriting 
( |3.9| ) . A manifestation of this ambiguity is seen in ( |4.7] ) , where the two extra unconstrained 
parameters in A" represent the flat directions of the scalar sector. On the other hand, the 
relation between % and the physical charges must be the same throughout the flow, as 



follows from ( 4.12| ), so that an attractor analysis would be sufficient for this purpose. In 



any case, one can always dualise the stabilisation equations for the seed solutions above to 
find any other solution and we comment on a possible way to construct % at the end of 
the next section. 

5. The static limit 

In this section we specialise the results of section ^ to the static case, using the ansatz of 
section and connect to the fake superpotential formalism. 

5.1 The static flow equations 

The static limit of the results in section ^ leads to several simplifications, since the solutions 
are necessarily spherically symmetric. This implies that w = and all quantities depend 
only on the radial coordinate. Similarly as for the actual field strength, spherical symmetry 
implies that the modified field strength IF is of the form: 

P = sine(ie ^dip®f , (5.1) 

where now F G H'^{Mcy) is fibred along r. By (3.9), it must reproduce the same black 
hole potential T^h as the physical charge F: 

l\\rf = Vb^ = l\\rf. (5.2) 

In this setting one chooses the arbitrary function e'" to be the phase of (F, fi). In terms 
of the inverse radial coordinate r = 1/r the first-order equations ( |3.7| , \i.8\) with r] = rjdr 
reduce to^ 

2 dr Im(e-^e-'"J7) - 2 Re(e-^e-'"0) = -f , (5.3) 
rj = -a-Qr = -Im{{t,n)/Z{r)), (5.4) 



^The signs depend on the conventions chosen for the Hodge dual. 
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where the second relation turns out to be equivalent to the second constraint ( |3.10| ). Ob- 
serve again that the presence of a nontrivial r] is essential for the generalisation of Denef's 
formalism with a fake field strength that is not a closed form. 



As described in section 3.2, at least in principle these equations can be further sim- 
plified by eliminating fake charges T from them and by using an ansatz J' for the scalars, 
whereupon we obtain equations for J'. In particular, since we already have an ansatz 



( |4.12 ) for the stabilisation equations, we can determine U and from it, so that eq. 



becomes: 

2 dr Im(e-^e-'"0) = drJ , (5.5) 
-t = drJ -2r]Re{e-^^e-'''nj), (5.6) 

where rj is given by the static limit of ( 4.24| ) as 



T] 



(?^^ {drTl,n) . (5.7) 



Recall from the previous section that our ansatz is automatically a solution of the constraint 
( 3.10D (and ( |5.4| )), and the equations of motion are solved if one can find a J , along the 



lines of section |4.2| , such that T constructed above reproduces the black hole potential in 
( |5.2| ), which now reads: 

\\\drJ\? = \\\Tf + (?^^{n,drnf . (5.8) 

This quadratic constraint can be used to relate the physical charges to the harmonic func- 
tions %, in addition to the generic requirements of section |4.2| . In summary, the static 
equations of motion are integrable if there exists an T-L and its corresponding TZ, constructed 



along the lines of section O, satisfying ( ^.8| ) . 



This is similar in spirit, but different than the approach of [|41| , f42| , |43| ], were one seeks 
to rewrite the black hole potential in (|5.2[ ) through a function of the physical charges 
and moduli z"" directly. In contrast, (^]^) is an equation relating the harmonic functions 
controlling the physical charges to the ones controlling the scalars through the period vector 

n. 

We have checked that J for the known explicit static solutions are such that they 
satisfy (5^) and hence are described by the flow equation (5^) with rj as in ( 4.24| ). Since 



all static non-BPS solutions are related by symplectic rotations to the seed solutions of 
section it follows that they satisfy the same duality-covariant equations. The nontrivial 
r/ is reflected in the anharmonic part of ( [4. 3D , controlled by the constant h that remains 



after setting the angular momentum to zero in ( ^4.11) . This observation is in line with |17|, 
where it was stressed that the crucial departure of the static non-BPS seed solution from a 
BPS-like ansatz is the presence of a parameter related to the asymptotic scalars, identified 
with this residual constant. 

5.2 The fake superpotential 

One can adopt the opposite point of view to the one taken above: first find fake charges T 
reproducing the black hole potential (^]^) and then solve the differential equations. Taking 
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the intersection product of both sides of ( ^.S] ) with the basis elements then leads to the 
equations for the scale factor and the scalars, which, with appropriate identifications, have 
the form of non-supersymmetric flow equations generated by a superpotential W p8| , 
analogous to the supersymmetric flow equations governed by the absolute value of the 
central charge: 

il = -e^e-'°(f , n) = -e^W , (5.9) 
i"^ = -e^e-'"5'^^(f, Vj,n) = -2e^g''^d-^W . (5.10) 

Whenever W is explicitly known for a given model and charge configuration, a practical 
way to connect it with our approach may be to first look for a moduli-independent matrix 
S that rotates the usual charge vector F so that: 

\z\:=\{sr,n)\ = w . (5.11) 

Then, its relation with F is defined by: 

f := iZn - ig^'^DaZVbn + ig^^VaZV-^n - iZCl , (5.12) 

where "DaZ = daZ + ^daKZ, with K being the Kahler potential. Note that in general 
F 7^ 5F, if we allow for S to be complex, and in fact this turns out to be the simplest 
choice. 

One can find the matrix 5* explicitly for the electric configuration, as in ( ^.3[ ) , assuming 
all physical scalars to have the same phase, say /. The relevant superpotential was given 
in |l^. Then, a suitable matrix S satisfying ( ^.11[ ) and defining F through ( 5.12| ) is 

S = diag(e"2'^,l,l,l,e2'^, 1,1,1) . (5.13) 

In terms of the parameters appearing in the solution of section ^ one can identify cot / = 
e^^M and check that the equations of motion ( [5. 3D are satisfied. The one- form rj is given 
by O). 

In the non-supersymmetric axion-free case M vanishes and / = 7r/2, so that S is 
constant (but not identity), while allowing for a r-dependent / leads to more general 
non-supersymmetric solutions. It is worth noting that r] = whenever S is constant 
(cf. eq. (5^)). In particular, when S = 1 we recover the supersymmetric case. 

Alternatively, one can rewrite (|5.6D and (|5.8| ) in terms of a real matrix T such that: 

TF := Ft = f - 2r] Re(e-^e-'"0) , (5.14) 

i(FT,oFT) = ybh + e-2V, (5.15) 
e-^"(FT,Jl) = VF-ie-^r/, (5.16) 

where W is the superpotential in ( |5.9D -( ^riO| ) . If T is known, it leads to simpler equations 
of motion for the scalars, that is 

2dr Im(e-^e-'"!^) = -Ft , (5.17) 



^°By the shorthand ST we mean rotating the symplectic vector of charges corresponding to the coefficients 
of F: S-ip',qjf, and arranging the resuh again as a three-form. 
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which have the advantage of being directly integrable to ( 3.24 ), giving the sought solutions 
as explained above. For the electric example above, T takes the form: 




c-^^ cot / 2qa" 

W= I ^ I , (5.18) 

^ 



so that (4^) can be written as J = —TTic, if the constants in (|4.2| ) are appropriately 
chosen. Similarly to its complex counterpart S, it reduces to a constant matrix in the 
axion-free case. 



It is interesting to point out that the matrix ( ^.18 ) is a (spacetime-dependent) element 



of the Peccei-Quinn group of transformations, defined as the largest subgroup of the sym- 
plectic group leaving the X^'s and the Kahler potential invariant. As was shown recently 
j66| j, applying such a transformation on the charges indeed shifts the black hole potential, 
as in ( 5.15[ ). For generic charges and phases of the scalars, the corresponding T can be 



found from the one in ( |5.18| ) by conjugation with the appropriate element of the symplec- 
tic group. Such a matrix would leave a certain combination of X^'s and Fj^s unchanged, 
e.g. for the magnetic dual of the electric solution in ( [4. 21) it would leave the F/'s invariant. 
Identifying the combinations that must be invariant for a given set of charges could be a 
way to determine T from first principles. 



6. Conclusions and outlook 

In this work, we have extended the formalism of |^, deriving symplectically covariant 
flow equations for non-BPS extremal black holes in = 2 supergravity, and we have 
constructed an ansatz for the corresponding stabilisation equations. The main novelty was 
to rewrite the electromagnetic part of the action in terms of a 'fake' field strength two- 
form ^ that does not have to be closed, where the non-closure turns out to be governed 
by a single one- form ij. The presence of this one- form is further related to the axions in 
the full black hole solution, apparently rendering such a deformation essential in a general 
description of non-BPS extremal black holes. Unfortunately, this complication makes the 
full equations challenging to solve directly, at least without considerable intuition about 
the form of the solution. 

To obtain that insight, we considered the known seed solution for under-rotating ex- 
tremal black holes in theories with cubic prepotentials. We showed that it can indeed be 
written in terms of stabilisation equations for the period vector, just as BPS black holes. 
The crucial difference is that the scalars are not stabilised in terms of harmonic functions 
only, but one finds that a ratio of harmonic functions is required. When the angular mo- 
mentum is set to zero, one simply has the inverse of a harmonic function, which vanishes 
near the horizon, but mixes with the other asymptotic constants at infinity. 

As one might expect, a comparison of these explicit solutions with our flow equations 
reveals the non-closure of JF^ to be reflected in exactly these non-harmonic parts of the 
stabilisation equations. Based on this, we proposed an ansatz for the generic stabilisation 
equations of under-rotating extremal black holes, satisfying several requirements coming 
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both from general arguments and known explicit solutions. Its practical realisation depends 
heavily on the model and in particular on invariants constructed from two charge vectors, 
one of which must correspond to a small black hole. Such invariants have been considered 
recently in |67, 63 1 in the context of multi-centre solutions. 

In the static case, we showed how to combine this ansatz with the first-order flow 
equations to identify the structure of J- for known solutions and infer its general form. It 
turned out that there are two ways of connecting the result to previous work. One involves 
a complex matrix resembling the matrix S introduced in the superpotential formalism |^8| . 
The other formulation simplifies the equations of motion, using the real matrix in ( 5.1 
belonging to the group of Peccei-Quinn transformations. This hints towards the possibility 
of obtaining such matrices systematically. 

Irrespectively of the precise description, one can characterise any solution by a vector 
of harmonic functions such that the quartic invariant computed on their poles is related 
to the one associated to the physical charges by a sign flip. In principle, it is possible to 
construct this vector using purely algebraic methods, which is equivalent to solving the 
attractor equations for static non-supersymmetric black holes. Once this is known for 
particular charge vectors, one can replace the vector of charges with harmonic functions to 
find 'H in (l4.12| ) and directly construct new solutions using our first-order equations. 

We expect this to generally hold also for the single-centre rotating solutions covered 
by our ansatz, as the only difference with respect to static solutions is in the choice for 
the angular momentum harmonic function, without modifying the structure of ( 4.12| ). 
In line with this expectation, we have verified that the proposed ansatz does satisfy the 
constraint (3.31). However, it is more difficult to impose ( p.9D , so it would be useful to find 
a generalisation of the arguments in section ^ and/or an extension of the fake superpotential 
formalism to the rotating case. 

It is important to note that our flow equations are by construction fully covariant with 
respect to electric-magnetic duality and are compatible with the general seed solutions 
in four dimensions. It then follows that they capture the full orbit of non-BPS extremal 
solutions for suitable choices of T ^ regardless of the existence of other stationary points 
of the action, which should not be part of the standard non-BPS orbit of extremal black 
holes. It is also interesting to point out the similarity to special cases of static non-extremal 
black hole solutions, which can be obtained through a deformation of extremal solutions 
controlled by a ratio of harmonic functions [38|, except that there it appears in the line 
element. 

On the microscopic side, it would be very interesting to reproduce the stabilisation 
equations ( [4.12| ). In the rotating case the ratio of harmonic functions survives the near- 
horizon limit and modifies the attractor equations, similarly to [22|, so one generally ex- 
pects this structure to be accessible from microscopies. Given the model of |^, where the 
constant part of M in (4.4) is interpreted as the angle between wrapped D3 branes, one ex- 
pects that the full angular momentum harmonic function might have a similar microscopic 
analogue. 

Finally, it is worthwhile stressing that albeit the explicit solutions that we have dis- 
cussed are only single centre, we have not made any assumptions on the number of centres 
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in the derivation of flow equations in section |^. Also the ansatz ( 4.12| ) is compatible with 
multi-centre harmonic functions. It would be illuminating to make a detailed compari- 
son with the results of |27, 28, 3C], as a test on the robustness of the assumption on the 
existence of stabilisation equations for generic extremal backgrounds. 
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A. On inverse harmonic functions 



It is possible to give a generic heuristic argument for the presence of ratios of harmonic 
functions of the kind seen in (O) in the solution for the scalars. We find it convenient to 



choose the arbitrary function e'° according to (|3.8|) , so that (|3.14D reduces to 

c = {g,g). (A.i) 

Similarly to the gauge part of ( |3.3D , this can be interpreted as an action for the tensor 

^ = J" - 2 Im ★oDfi + 2 Re D(e^^J7) A {dt + oj) , (A.2) 
which respects the same pseudo-selfduality condition (2.14) as J^, assuming that T does: 

Mi<g = g. (A.3) 
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The scalar part can be shown to be pseudo-selfdual using (|2l3|) and The matrix M 



is crucial for the existence of such a constraint, since it is not possible to impose selfduality 
on a four-dimensional field strength unless one complexifies it. However, it can be done in 
4n + 2 dimensions, and ( [A.3| ) descends from the ten-dimensional constraint on the five- form. 

For the ordinary Einstein-Maxwell theory, gauge field equations in backgrounds of the 
type ( |3.lD naturally lead to harmonic functions (cf. eg. appendix B in |61]). Motivated by 



( |4.5| ), we further assume that the rotational one-form satisfies duj = ★od-^^- Denote spatial 
directions by i = 1, 2, 3 and consider first electric solutions, Fij = 0, for which the Bianchi 
identity implies the existence of an electrostatic potential 

M 

diFtj = djFu Fu = d,— . (A.4) 

The field equations further impose that H is harmonic 



e 



-^ = H, V^H = 0. (A.5) 



Similarly, for a magnetic solution, Fu = 0, the Bianchi identities relate the field strength 
to a harmonic function 

e''^''diFjk = ^ Fij = e,,kdkH, (A.6) 

where e^k is the Levi-Civita permutation symbol and the Einstein equation implies again 
= H. These solutions are related by an electric-magnetic duality rotation and belong 
to a class of solutions called the Majumdar-Papapetrou solutions. 

Now, consider the case that the field strength is constrained to be pseudo-selfdual, 
as in ( |A.3| ). Then, the distinction between equations of motion and Bianchi identities 
disappears, and the solutions can no longer be purely electric or purely magnetic. For such 
a field, the gauge part of the action in ( |3.3| ) leads to the equation of motion (see [53| for 
details) 

d(^-7W*J-) = 0, (A.7) 

which is then solved for both cases: 

J" = -ko^H , e"^ = I{H) , (A.8) 

^ = *„df, = (A.9) 

where now H denotes a vector of harmonic functions and / is a model-dependent invariant. 
Imposing closure of which is equivalent to the existence of a gauge potential in four 
dimensions, one concludes that only the first solution survives, leading to the standard 
description by harmonic functions. 

Going back to the Lagrangian in (|A.l| ), the above discussion of selfdual fields applies 
including the scalar sector. In view of this, unlike IF, is not necessarily closed and the 
two independent solutions in (|A.8D -( |Al9| ) are allowed. Based on this, one concludes that 
a vector of inverse harmonic functions is a zero mode of the equations of motion following 
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from ( [A.l| ). Such a vector must be part of the general solution for the scalar sector only, 
given that the gauge field T is described by harmonic functions, as above. 

Since the equations of motion are nonlinear, it is a nontrivial task to find full solutions 
with both kinds of zero modes turned on. Nevertheless, this discussion demonstrates that 
one can consider the scalar sector a priori independently from the vector one. It also lends 
credibility to the presence of ratios of harmonic functions in the ansatz ( [4.12 ) and the 
different sets of harmonic functions in (4.3). 
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